Metal-Insulator Transition Accompanied with a Charge Ordering in the 

One-dimensional t — J' Model 
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We study the metahinsulator transition accompanied with a charge ordering in the one-dimensional 
(ID) t — J' model at quarter filling by the density matrix renormalization group method. In 
this model the nearest-neighbor hopping energy t competes with the next-nearest-neighbor ex- 
change energy J'. We have found that a metal-insulator transition occurs at a finite value of t/J'; 
{t/J')c — 0.18 and the transition is of first order. In the insulating phase for small t/J' , there is an 
alternating charge ordering and the system behaves as a ID quantum Heisenberg antiferromagnet. 
The metallic side belongs to the universality class of the Tomonaga-Luttinger liquids. The quantum 
phase transition is an example of melting of the ID quantum Heisenberg antiferromagnet. 
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Recently low-dimensional copper oxides have attracted 
much attention as a route to understand the physics of 
the high-Tc superconductivity. Since the superconduc- 
tivity in a doped ladder system with a spin gap was 
theoretically suggested ||], several experimental stud- 
ies of the quasi one-dimensional (ID) systems with the 
ladder structure have been carried out. Especially, the 
discovery of superconductivity in the quasi ID material 
Sra;Cai4_a;Cu2404i.84 {x — 0.4) undcr a high pressure 
has received considerable interest The parent mate- 
rial Sri4Cu2404i includes Cu02-chain layers and CU2O3- 
ladder layers separated by the Sr layers Spin exci- 
tation gaps are observed in both chains and ladders in 
ambient pressure 

[ill- The origin of the spin gap in 
the chain is still controversial in spite of intensive studies 
|^-|ll|. By neutron scattering measurements, Matsuda et 
al. concluded that the spin gap in the chain originates 
from dimer formation between the copper pairs separated 
by two and four distances of the nearest neighbor cop- 
per ions. On the other hand, in more recent neutron 
scattering measurements Eccleston et al. suggested 
that dimers are formed between next-nearest-neighbor 
Cu ions which are arranged with a period of 5 Cu-Cu dis- 
tances in the chain . In the both models considered by 
Eccleston et al. and Matsuda et al., the superexchange 
interactions between the next-nearest-neighbor Cu spins 
play an important role. 

Matsuda et oZ.'s results have been interpreted in ||] 
that there exists the alternating Cu"*"^ and Cu"*""^ order 
along the chain. According to this scenario, the mean 
valence of Cu ion in the chain is Cu"*"^'^. On the con- 
trary, Cu"^^'^ is assumed in Eccleston et aZ.'s results 
Which scenario is realized is still an open problem. Note 
that for either case of Cu+^-^ or Cu+^ '', the Cu 3d band 
is partially filled; the highest 3d band is nearly quarter- 
filled. Although the Cu 3d band is partially filled, the 
electrical resistivity of the system shows a semiconduct- 
ing temperature dependence [H. Precise understanding 



of the semiconducting temperature dependence depends 
on whether the system is just at the quarter filling or not, 
since localization effects should be taken into account for 
the latter case. 

In the edge-sharing Cu02 chain, the superexchange 
interaction between next-nearest-neighbor Cu spins 
through Cu-O-O-Cu path is expected to be more im- 
portant than the nearest-neighbor superexchange inter- 
action, since in the latter case the bond angle of Cu-O-Cu 
is close to 90 degree. Thus, the model we consider in the 
present study is the simplest case where the next-nearest- 
neighbor superexchange interaction competes with the 
kinetic energy. For the kinetic energy we consider only 
the nearest-neighbor hoppings, since in the vicinity of the 
metal-insulator transition of a quarter filled system the 
nearest-neighbor hoppings are more important than the 
next-nearest-neighbor hoppings even if the matrix ele- 
ment of the latter processes is bigger than the former. 
We assume strong Coulomb repulsion between the d- 
electrons, leading to suppression of the double occupancy 
of the carriers. This model may be called as the one- 
dimensional t — J' model where t represents the nearest- 
neighbor hoppings and J' the next-nearest-neighbor ex- 
change. 




FIG. 1. The t — J' chain model. Each dot denotes a site 
and arrows represent carriers with spins specified by arrows. 



The model Hamiltonian is defined as 

n = ^(4^+1, + H.C.) + J'Y,S^■ S,+2, (1) 
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where the operator c^^. (cio-) creates (annihilates) a parti- 
cle with a spin a at ith site. These operators are defined 
in the subspace where double occupancy is excluded. Si 
denotes the spin operator at ith site. The model is topo- 
logically equivalent to the zigzag chain which is shown 
in Fig. The solid and broken lines denote the nearest- 
neighbor hoppings t and the next-nearest-neighbor ex- 
change interactions J', respectively. 

In this paper we restrict ourselves to the quarter fill- 
ing. At this filling, the system has two simple limits, (i) 
t 0, J' is finite : In this limit, all particles localize 
at odd or even sites, and the system has a gap in the 
charge excitations. This charge-ordered state behaves as 
an antiferromagnetic (AF) Heisenberg chain with the ex- 
change interaction J' ||lj]. In the AF Heisenberg chain, it 
is well known that there is no gap in the spin excitations, 
(ii) J' ^ 0, t is finite : In the case of J' — 0, the system 
is identical to the spinless fermions with the macroscopic 
spin degeneracy of 2^ {N is the number of spins). For an 
infinitesimal J' the spin degeneracy is lifted and the sys- 
tem is expected to belong to the universality class of the 
Tomonaga-Luttinger (TL) liquids In this limit, the 
system is gapless in both the charge and spin excitations. 
Since the present system has the two limits mentioned 
above, one can expect that a quantum phase transition 
occurs at some critical value {t/J')c from an insulating 
phase (the charge gap Ac is finite) to a metallic phase 
(Ac = 0) as increasing t/J' from to oo. Since parti- 
cles are expected to order alternatively for a smaller t/J' 
than (t/J')c, the quantum phase transition is a charge 
ordering transition. 

The purpose of the present study is to answer the fol- 
lowing two questions: (i) Does the quantum phase tran- 
sition occur at a finite (t/ J')c or already at infinitesimal 
t/J' 7 (ii) What is the nature of the phase transition? Is 
it a first- or second-order transition? 

In order to determine the transition point, we calcu- 
late two quantities: the charge density order parame- 
ter ricd and the charge gap Ac. In the insulating phase, 
it is expected that the all particles localize at even or 
odd sites. Therefore, we consider as the order parameter 
n-cd — "odd — "■oven, whcre nodd(f^cvon) is the mean value 
of the number density rii at odd (even) sites. The order 
parameter ricd is defined as 

' i=odd 

While ricd should be zero in the metallic phase, it should 
be finite in the insulating phase. In a finite-size system, 
the symmetry breaking does not occur if there is no ex- 
ternal field. Thus, we calculate ricd with a small external 
field e at the end of the system. This external field is 
added to the model Hamiltonian (||) as H — ej^a Ci<jC1(t 
at the first site (e > 0). We study the response of the 
system against the small external field. 
The charge gap is defined by 

2Ac(i) = EgiNo + 2, 0; L) + E^iNo - 2, 0; L) 



-2Eg{No,0;L), (3) 

where Eg{N, Sz', L) denotes the ground state energy of 
the system in which the number of particles, the z com- 
ponent of the total spin and the system size are denoted 
by N, Sz and L, respectively. In the quarter-filled case, 
N = No = L/2. 
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FIG. 2. Local number density Ui at each site for various 
values of t/J'- t/J' = 0.05, 0.10, 0.13, 0.14, 0.15 and 0.17. 
For a small t/J', riodd — 1 and ncvcn — 0. Ths system size L 
of this example is 32 (e/J' = 1.0 x 10"^). 

To study the quantum phase transition numerically, 
it is essential to look at the size dependence system- 
atically and therefore we need to treat large-size sys- 
tems. To this end, we use the density matrix renormal- 
ization group (DMRG) method jl^], which is one of the 
standard numerical methods to study one-dimensional 
quantum systems. We study systems of various sizes 
(L = 12, 16, 24, 32, and 48) by using the finite-system al- 
gorithm of the DMRG method ||l|l. We keep 100 - 150 
states in the renormalization procedures so that trunca- 
tion errors in the density matrix are less than 1.0 x lO^"' 
(typical numerical errors are about 1.0 x 10^^). 




o; 1^^^ : 
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FIG. 3. Results of size extrapolation of n^d vs t/J' (see 
the text). In the inset we plot n^d for various system sizes. 

Figure ^ shows examples of the local number density 
Hi for various values oit/ J' with an external field e/ J' = 
1.0 X lO"'^ (L — 32). Shown here are central 28 sites 
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excluding 4 sites at both ends of the system. For small 
t/ J'{<Si 1), one can clearly see charge ordering behaviors. 
Most particles spontaneously localize at odd or even sites 
and in this example the odd sites are selected by the small 
symmetry breaking field at the first site. For a small t/ J' 
charge fluctuations are small, i.e., riodd — 1 and nevcn — 
0. As increasing the value of t/J', the difference between 
JT-odd and ricvcn decreases for a fixed system size. Both 
jT-odd and novcn approach the mean value of the number 
density N^/L = 0.5 for large t/ J' . 

Charge density order parameters n^d are shown in the 
inset of Fig. || for various system sizes {L = 16,24,32, 
and 48). Error bars in the inset represent the standard 
deviation. 




cd' 



(4) 



but for most of data they are smaller than the symbols. 
One can see that as the system size is increased, the 
disappearance of ricd becomes sharper. The region of the 
charge ordered phase spreads out for larger-size systems. 
On the other hand, for all system sizes Tied vanishes in 
the region t/ J' > 0.2. These results may be an indication 
that the charge ordering transition is of the first order. 
However to draw a definite conclusion careful finite size 
scaling is necessary. Results of size extrapolation of ricd 
are shown in the main part of Fig. |^. We have carried 
out the linear extrapolation using data for the system 
size L = 24 and 32. For the transition region 0.13 < 
t/ J' < 0.20, we have used data for L = 32 and 48. In the 
case that the extrapolated value obtained by the simple 
linear extrapolation is greater than 1, the upper error 
bar is fixed to unity. The values of Tied for t/ J' ^ 0.2 are 
extrapolated to zero, while the values below t/J' ~ 0.2 
tend to be finite. Concerning the symmetry breaking 
field, we have confirmed that ricd in the metallic side 
are suppressed even more and the disappearance of rtcd 
becomes sharper for a smaller external field e (= 1.0 x 

10-V). 
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FIG. 4. Size dependence of the cliarge gap Ac {L) scaled by 
t for several values of t/J'; t/J' = 0.14, 0.16, 0.18, 0.20, 0.22, 
0.24 and 0.26. Note that the system size dependence is qual- 
itatively different between t/J' = 0.14,0.16 and t/J' > 0.18. 



Now we turn to the charge gap. In Fig. |[ we present 
the system-size dependence of the charge gap scaled by 
t for several values of t/J'. The system-size dependence 
of Ac{L) for t/J' < 0.16 is clearly different from that 
for t/J' > 0.18. For t/J' < 0.16, Ac(i) increases as L 
is increased, which suggests that the infinite-size limit of 
Ac is finite. One can conclude that the system is in the 
insulating phase for t/J' < 0.16. On the other hand, for 
t/J' > 0.18, it tends to zero in the infinite-size limit; the 
system is in the metallic phase. The observed behaviors 
of the charge gap and the system-size dependence of the 
charge density order parameters are compatible with the 
first-order phase transition at {t/J')c ~ 0.18. 
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FIG. 5. Size dependence of the spin gap As (I/) scaled by 
J' for several values of t/J'; t/J' = 0.14, 0.16, 0.18, 0.20, 
0.22, 0.24 and 0.26 (from top to bottom). The inset shows Xs 
evaluated from As{L) (see the text). 



We also calculate the spin gap defined as 
As(L) = E^iNo, 1; L) ~ E^{No,0; L), 



(5) 



in order to see whether there is an anomaly in the 
spin excitations at the quantum phase transition. The 
system-size dependence of the spin gap scaled by J' 
is shown in Fig. ^ In both (insulating and metallic) 
phases, As(L) tends to zero in the infinite-size limit and 
one can conclude that there is no gap in the spin ex- 
citations. The slope dAs(L)/d(l/L)|(i/i)=o gives the 
spin susceptibility Xs of the system through the relation 
As(L) = 2/{xsL) [|l6[. The coefficient of the linear term 
is determined by a fitting to a quadratic function includ- 
ing the origin;(l/L = 0, As{L — oo) — 0), points for 
L = 48; (l/48,As(48)) and for L = 32; (l/32,As(32)). 
We obtain the coefficients also by using cubic and bi- 
quadratic functions including four (the origin, points for 
L — 48, 32, and 24) and five (the origin, points for 
L — 48, 32, 24, and 16) points, respectively. The inset of 
Fig. H shows Xs ys t/J' obtained by the above-mentioned 
procedures. Good coincidence between the results evalu- 
ated from cubic and biquadratic functions indicates that 
the values thus obtained for Xs are reliable. One can 
see that t/J'-dependence of Xs for t/J' < 0.18 may be 
slightly different from that for t/J' > 0.18, although it 
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is not possible to conclude whether there is a discontinu- 
ity at (t/J')c — 0.18 or not. Since the phase transition 
occurs at J' which is about five times bigger than t, it 
is natural to assume that the energy scale of system in 
the region near the critical point is dominated by J'. A 
natural consequence is that the anomaly at the critical 
point, if any, should be small. 

The present model is simple but rather special in 
the sense that the nearest-neighbor exchange energy J 
and the next-nearest-neighbor hopping energy t' are ne- 
glected when t and J' are considered. In spite of the spe- 
cial nature of the model, the metal-insulator transition 
with the charge ordering is generic. To show this, we in- 
vestigate a few cases with the nearest-neighbor exchange 
or the next-nearest-neighbor hopping. In Fig. ^, we show 
Ucd in the following two cases: with an additional nearest- 
neighbor (ferromagnetic) exchange J/ J' = —0.5 and with 
a next-nearest-neighbor hopping t' — t. In both cases, 
for small t/ J', ^cd is nearly equal to 1, and it disappears 
rapidly around some critical value oi t/ J'. These behav- 
iors are qualitatively similar to that of the t — J' chain 
model which we have discussed in the present paper in 
detail. In the case with J [t — J — J' model) , the criti- 
cal value of t/ J' shifts to a larger value than that of the 
t — J' chain model since the hopping energy t competes 
with not only J' but also J. On the other hand, in the 
case with t'{— t) {t — t' ~ J' model), the critical value of 
t/ J' is smaller than that of the t — J' chain model. This 
is natural since in the t — t' — J' model the gain of the 
kinetic energy in the metallic phase is larger than that in 
the t — J' model. 
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FIG. 6. Comparison of n^d in the t — J — J' chain 
(J/ J' = -0.5) and t ~ t' - J' chain {t' = t) with that in the 
t- J' chain. The system size is L = 32 (e/J' = 1.0 x 10"^). 

In conclusion, we have calculated the charge density or- 
der parameter, the charge gap, and the spin gap obtained 
in the t — J' chain model using the DMRG method. We 
have found that the quantum phase transition from a 
charge ordered insulator to a metal occurs at t/J' ~ 0.18 
at quarter filling, and the transition is of first order. The 
transition is generic in the models in which the kinetic en- 
ergy competes with the next-nearest-neighbor exchange 
energy at quarter filling. For the present model, an in- 
teresting question is existence or absence of the spin gap 



close to the metal-insulator transition. The spin excita- 
tions are gapless in both limits t — > and J' ^ 0. We 
have found that the spin excitations remain gapless over 
the whole t/J' range between the metallic side and the 
insulating side. The spin susceptibility is obtained ^from 
the spin gap of the finite systems. One can see that the 
t/ J'-dependence of Xs in the metallic phase is slightly 
different from that in the insulating phase, and there is 
a possibility that a small anomaly exists at the critical 
point in the bulk system. 

The quantum phase transition discussed in the present 
paper is an example of melting of the ID quantum 
Heisenberg antiferromagnet. In this example, the melt- 
ing is of first order and the metallic side belongs to the 
universality class of the TL liquids. 
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